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We analyze the dynamic toroidal multipoles and prove that they do not have an independent 
physical meaning with respect to their interaction with electromagnetic waves. We analytically show 
how the split into electric and toroidal parts causes the appearance of non-radiative components in 
each of the two parts. These non-radiative components, which cancel each other when both parts 
are summed, preclude the separate determination of each part by means of measurements of the 
radiation from the source or of its coupling to external electromagnetic waves. In other words, there 
is no toroidal radiation or independent toroidal electromagnetic coupling. The formal meaning of 
the toroidal multipoles is clear in our derivations. They are the higher order terms of an expansion 
of the multipolar coefficients of electric parity with respect to the electromagnetic size of the source. 


I. INTRODUCTION, SUMMARY AND 
OUTLINE 

After the introduction of toroidal multipoles in the 
static case [lj, the dynamic toroidal multipoles were pre¬ 
sented as a new independent multipole family that had 
been previously ignored ai- The new family was de¬ 
rived from the split of the transverse multipoles of elec¬ 
tric parity into two parts. These parts are often referred 
to as electric and toroidal. The dynamic toroidal multi¬ 
poles have also been analyzed in Refs. 0 and a and are 
currently being considered in the areas of metamaterials, 
plasmonics, and nanophotonics am. 

In this article, we analyze the split between the dy¬ 
namic electric and toroidal multipoles of a localized 
source distribution. We find that the two parts cannot be 
separately determined by measuring the electromagnetic 
fields produced by the source outside the source region. 
The two parts can also not be separately determined by 
measuring the coupling between the source and exter¬ 
nally incident electromagnetic waves. The result applies 
to both near and far field situations, and implies that 
there is no independent coupling of toroidal character 
between the sources and the electromagnetic field, and 
that there is no radiation of pure toroidal character. The 
toroidal multipoles can hence not be considered an in¬ 
dependent family. Rather, our analysis makes clear that 
the electric and toroidal parts correspond, respectively, 
to the lowest and the higher order terms of the expansion 
of the exact multipolar coefficients of electric parity with 
respect to the electromagnetic size of the source. In this 
respect, our article contains analytical expressions that 
inherently contain all correction orders: They are exact 
for any source size. 

Our results do not question the usefulness of consid¬ 
ering the higher order terms. For example, some ex¬ 
perimental results can only be explained by adding the 
toroidal dipole contributions to the lowest order terms 
(see Fig. 3c in Ref. 0, and Fig. 4 in Ref. HI)- The 
same is true for non-radiating spherical configurations 
(see Fig. 2 in Ref. [20h . These cases are examples where 


the toroidal terms are the dominant dipolar terms. 

In the following, we first outline the mathematical 
setting in which we carry out the analysis. We then 
highlight the direct connections between our momentum 
space approach and the physics of multipolar couplings 
between sources and fields. On the one hand, three kinds 
of degrees of freedom are needed to describe the sources. 
Two of them are transverse and have opposite (electric 
vs. magnetic) parity, and the other one is longitudinal. 
On the other hand, only the two transverse degrees of 
freedom are needed to describe the electromagnetic fields 
produced by the sources outside the source region, where 
the field produced by the longitudinal degrees of freedom 
is identically zero. 

Then, we analyze the physical meaning of the toroidal 
multipoles. Our analysis starts from recently obtained 
exact analytical expressions for the transverse multipo¬ 
lar coefficients of electricparity. These expressions are 
valid for any source size J21| . We take the dipolar case as 
the guiding example, make the small source approxima¬ 
tion keeping the two lowest orders, and show that: The 
lowest order term is the well know approximate expres¬ 
sion for the electric dipole of electromagnetic.ally small 
sources; the term of second lowest order is the approx¬ 
imate expression for the toroidal dipole of electromag- 
netically small sources [Eq. (2.11) in Ref. 3}. We then 
analyze the original split [2], [3| in the general multipolar 
case, and establish that both electric and toroidal parts 
contain non-radiative components, which appear due to 
the separation of terms of different order in the size of 
the source. These non-radiative components, which can¬ 
cel each other when both parts are summed, preclude 
the separate determination of each part by measuring 
the radiation from the source outside the source region, 
or by measuring the coupling between the source and ex¬ 
ternally incident electromagnetic waves. Afterwards, we 
analyze a different class of splits where the two parts are 
free of this kind of non-radiative components. We find 
that, in this case, what precludes their independent mea¬ 
surement is the absence of longitudinal radiation. We 
show in App. [C] that this is the same reason that pre¬ 
cludes the independent measurement of the recently pro- 
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posed re-definition of toroidal multipoles (Box 2 in Ref. 


Finally, we use the obtained insights to clarify some 
statements that are often found in the literature. 


II. MATHEMATICAL SETTING 


We start with an electric current density distribution 
J(r,t) embedded in an isotropic homogeneous medium 
characterized by real valued permittivity e and perme¬ 
ability /i. We assume that J(r, f) is confined in space 
so that J(r, t) = 0 for |r| > R. We consider its energy- 
momentum Fourier representation 


J(r, t) = 1Z 
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and treat each frequency w separately. The lower limit of 
the integral in dui excludes the static case w = 0, which 
we do not treat in this paper. This formal setting is 
independent of the cause of the current. For example, 
in a scattering situation the current is induced in the 
scatterer by an external field. 



FIG. 1. The monochromatic electric current density distri¬ 
bution Ju(r) is confined to a sphere of radius R: The source 
region. In this article, we consider the electromagnetic fields 
produced by the source outside its region, i.e. [E„(r), H w (r)] 
for |r| > R. These fields are completely determined by 
the Fourier components of the current density which meet 
|p| — uo /c (see Fig. 0. The other components do not produce 
fields outside the source region. This is a general result that 
is also valid when the confining volume defining the source 
region is not spherical (see text). 


The electromagnetic fields radiated by the source at 
each frequency ui are solely determined by the part of 
J aj (p) in the domain that satisfies |p| = ujy/eJI = w/c. 
This theorem has been proven in Refs. [22] and The 
authors showed that this part of the current completely 
determines the electromagnetic field generated by the 
source at any point r such that |r| > R , which, in 
particular, includes near fields. This result has been 
extended beyond spheres to confining volumes of any 
smooth enough shape (see Chap. 9 in Ref. [24]): The part 
of J w (p) in the domain that satisfies |p| = Wy/ejl = w/c 
determines the electromagnetic field generated by the 
source at any point outside the source region , for a wide 
class of shapes of such source region. Figure [T| depicts 
a spatially confined monochromatic source distribution 
and the fields that it generates outside the source region. 

We denote by J w (p) the components of J w (p) in the 
spherical shell of radius |p| = w/c as a function of the 
solid angle p = p/|p| [see region (al) in Fig. 0. 

In this domain, the relevant scalar product is (A\B) = 
f dp A^(p)B(p), where t denotes conjugate transpose, 
and p runs over the entire spherical shell. 

We expand J u (p) in an orthonormal basis for vector 
functions defined in a spherical shell: 

■Up) = y " a jm^jm(p) T bj m X.j rn { p) + Cj m Wj Tn ( p). 

jm 

( 2 ) 

The basis is composed by the three families of vector 
multipolar functions in momentum space (Sec. B/.3 in 
Ref. [25|) 


Xjm(p) — r—-— Llj m (p), 

VJU +1) 

Zjm (p) — ^P X Xy m (p}, 

w im (p) = pY jm (p), 


( 3 ) 


where the Yj m (p) are the scalar spherical harmonics and 
the three components of the vector L are the angular 
momentum operators for scalar functions. 

Each of the functions in Eq. (0) is an eigenstate of the 
angular momentum squared J 2 , the angular momentum 
along one axis a, for which we choose a. = z, and the 
parity operator II. As depicted in Fig. 0 the polarization 
of X ?m (p) and Z jm (p) is orthogonal (transverse) to p, 
and the polarization of W jm (p) is parallel (longitudinal) 
to p [see region (a2) in Fig. 0. In coordinate space (r), 
this distinction corresponds to the distinction between 
divergence free (transverse) and curl free (longitudinal) 
vector fields. Table G] contains the eigenvalues for the 
aforementioned three operators and the transverse (_L) 
or longitudinal (||) character of each vector multipolar 
function. 

We provide expressions for the {a“ m , b^ m , c^ m } in App. 
E which we derived in Ref. I2II . 
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FIG. 2. Only the transverse components of the current den¬ 
sity J w (p) with p = produce electromagnetic fields out¬ 
side the source region (see Fig. [I}. (al) Spherical shell of 
radius |p| = — in momentum space. The (a2) region de¬ 
picts two different orthonormal bases for vector functions de¬ 
fined on the shell. The {X jm (P) 5 (p),W jm (p)} and the 
Y j,i+i,m(p),Yj,j_ 1 , m (p)}. The ( j,m ) indexes 
are omitted in the figure, and the subscript in the Y’s is the 
difference between their second and first indexes. The two 
bases are related to each other as written in Eq. (f7|l . With 
respect to their polarization vectors: X jm (P) 5 Zj7jx (p), and 
Yj,j,m(p) are orthogonal (transverse) to the momentum vec¬ 
tor p, Wj m (p) is parallel (longitudinal) to p, and Yj > 3 _|_i irrl (p) 
and Yjj_i, m (p) lay in the W jm (p) — Zj m (p) plane and are 
neither transverse nor longitudinal, as seen in (a3). The trans¬ 
verse vectors are represented by dashed blue arrows, the lon¬ 
gitudinal ones by solid black arrows and the ones of mixed 
character by dotted red arrows. For the figure, we took j = 1 
in Eq. Q for the relationships between (Z(p),W(p)} and 
{Y+i(p), Y_r(p)}. 
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TABLE I. Vector multipolar functions: Polarization character 
and eigenvalues of J 2 , J z and parity II. The symbol T means 
transverse polarization. The symbol || means longitudinal po¬ 
larization. Both j and m are integers, and m = — j .. .j. For 
X jm (p) and Z jrn (p), j > 0. For W jm (p), j > 0. 


A. Multipolar coupling between sources and fields 
in momentum space 

The decomposition in Eq. © is useful for discussing 
the physics of multipolar couplings between the sources 
and the fields. The complex scalars deter¬ 

mine the coupling of the source to the transverse elec¬ 
tromagnetic field, a.k.a radiation field, a.k.a electromag¬ 
netic waves. The one-to-one relationship between the 
coefficients of the source and the multipolar 
fields of frequency w radiated by it is exposed in App. B 
of Ref. |H, Chap. 9.10 in Ref. (13, and Ref. HD, for ex¬ 
ample, where the analysis is made in the r domain. The 
{dj m ,bj m } represent independent degrees of freedom in 
the following sense. Each individual coefficient is con¬ 
nected to the radiation (absorption) of a particular kind 
of electromagnetic wave: A transverse multipolar field 
with well defined frequency, angular momentum, and 
parity. The a^ m are connected to the multipolar waves 
of electric parity and the to those of magnetic parity 
(see Tab. [[]). The multipolar waves connected to differ¬ 
ent coefficients are orthogonal to each other. In prin¬ 
ciple, each individual {a" m ,fe“ TO } can be independently 
measured through the interaction of the source with the 
radiation field. 

The represent the longitudinal degrees of freedom 
of J^j (p) , but the longitudinal electric field with |p| = w/c 
is zero outside the source region defined in Fig. [T] This 
can be explained in the following way. The longitudinal 
components of the current (c" OT ) generate a longitudinal 
electric field outside the source region defined in Fig. |TJ 
but such field is canceled by the longitudinal electric field 
produced by the associated charge density. This state¬ 
ment follows from the continuity equation, as we show 
in App. [Bj The cancellation can also be recognized in 
§13.3 pl875-1877 in Ref. [28|, and in the expression of 
the electric field in Eq. (2.39) in Ref. [§:, where the two 
contributions to the longitudinal field cancel exactly. In 
Ref. [3 the longitudinal wave functions are ( r )- 

The bottom line is that the {a" m , 6" m } contain the nec¬ 
essary and sufficient information to determine both the 
fields of frequency ui produced by the sources outside the 
source region, and the coupling of the sources to external 
electromagnetic waves of such frequency. 

The lack of longitudinal radiation means that, even 
though there are three kinds of degrees of freedom in 
the current, {a“ m , i>" m , c“ m }, only two of them are ob¬ 
servable from outside the source region. The powerful 
theory developed by Gustav Mie for spherical scatterers 
[29j is probably the simplest example of this. In this 
case, the current density is induced in the sphere by an 
external incident field. The {a" m ,fr" m } of this induced 
current are determined by the Mie coefficients, which 
can be obtained analytically thanks to the symmetry of 
the problem. Using the Mie coefficients, the field pro¬ 
duced by the induced current can be determined exactly 
everywhere outside the sphere. This allows the exact 
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calculation of near fields, which is crucial in many appli¬ 
cations like for example: Determination of the transition 
rates of atoms placed near spheres [2Qj, determination of 
the forces and torques caused by evanescent fields onto a 
sphere near a substrate fU, modeling the tips of scan¬ 
ning tunnelling optical microscopes [321 ] , and general op¬ 
tical trapping analysis [33j ]. 

III. THE DYNAMIC TOROIDAL MULTIPOLES 

We will now study the split of the multipolar coeffi¬ 
cients of electric parity (a“ m ) which lead to the appear¬ 
ance of the dynamic toroidal multipoles 10- We use 
the dipolar case j = 1 as our guiding example and then 
generalize the findings to any j > 1. 

We first show that the well known expressions for the 
electric and toroidal dipoles of small sources are, respec¬ 
tively, the lowest and second lowest order terms in a series 
expansion of the exact expression of af m . 


approximation of the exact expression of a“. All higher 
order terms are easily obtained from the Taylor series of 
the spherical Bessel functions in Eq. 

The same structure revealed in Eqs. © to © under¬ 
lies the general case of the electric/toroidal split for j > 1. 
In Refs. @ and 0 the lowest order term in a" m , of order 
(fe) J_1 , is isolated as in Eq. ©. The rest, a" — e“ in our 
example, is called toroidal multipole form factor. Then, 
the lowest order term of the toroidal multipole form fac¬ 
tor, of order (fc) J+1 [like our Eq. ©], is called toroidal 
multipole moment. For simplicity, we refer to the toroidal 
multipole form factor as the toroidal part. It corresponds 
to the toroidal multipolar family as defined in Refs. l2jand 
0 In the dipolar example, the toroidal part (a“ — e“) 
is equal to —%tf plus all other higher order corrections 

7T V 3 

coming from the terms (fc) s>2 of the expansion of the 
spherical Bessel functions in Eq. ©• 

The formal origin of the toroidal multipoles is recog¬ 
nized in the literature 0, [l0 [I0 [2(| • We now investigate 
their physical meaning. 


A. Formal meaning 


We start from the recently obtained exact expression 
for the dipolar vector of electric parity in the spherical 
vector basis af = [a^, af 0 , af_f\ T (Eq. (21) in Ref. 


a 


UJ 

1 


d3r 

d 3 r {3 [f f J w (r)] f -J w (r)} j 2 {kr), 


(4) 


where k = u/c, and t means conjugate transpose. Let 
us now make the small argument approximation to the 
spherical Bessel functions with terms up to second order 
[ jo{kr) ps 1 — (fcr) 2 /6 and j 2 (kr) ps (kr) 2 / 15], and then 
group the contributions with the same power of k. 



e i 


~ fc2 J ^ { [r f J^(r)] r - 2r 2 J w (r)} . (6) 

^-v-' 

*1 

We find that the (fc)° term, e“ in Eq. ©. is the well 
known small source approximation of the electric dipole 
moment of a current density distribution [Eq. (9.14) in 
Ref. H3]. The k 2 term , t“ in Eq. ©. is the small 
source approximation of the toroidal dipole moment [Eq. 
(2.11) in Ref. 0. The integral in Eq. (6) contains 
contributions from the two integrals in Eq. (4): It is the 
sum of the terms of order k 2 coming from both jo(kr ) in 
the first integral of Eq. (4) and j 2 (kr) in the second. 

It is now clear that t“ is nothing but the next to lead¬ 
ing order term in the electromagnetically small source 


B. Physical meaning 

The splitting of a“ m into two parts 00 has a physi¬ 
cally relevant consequence. Both parts, often called elec¬ 
tric and toroidal, contain Fourier components of the cur¬ 
rent from outside the |p| = u/c shell. The |p| y^ 
components do not couple to electromagnetic waves 
For example, it follows from the properties of the Fourier 
transform that the result of the integral f d 3 rJ w (r) in 
the expression of the electric part in Eq. © is propor¬ 
tional to Joj(p = 0). Since |p| = 0 y^ w/c, it follows that 
the toroidal part, af — ef, must also contain a contribu¬ 
tion proportional to J w (p = 0). This is so because 
does not have any out of shell contribution [Eq. (4.12a) 
in Ref. [23| : The out of shell contributions in the elec¬ 
tric part must be canceled by those in the toroidal part 
when the two parts are summed. Since the |p| y^ ui/c 
components do not couple to the transverse electromag¬ 
netic field, it is impossible to determine the toroidal (elec¬ 
tric) parts of a“ by measuring the fields produced by the 
source outside its region, or by externally exciting the 
source with electromagnetic waves. As previously dis¬ 
cussed, this conclusion does not depend on whether the 
measurement (excitation) occurs in(from) the near, mid, 
or far field zones, as long as they occur outside the source 
region. For the sake of discussion, let us now consider 
the hypothetical case where the first term in the Taylor 
series of a given component of a“ is zero at a partic¬ 
ular frequency. We make no judgment about whether 
this hypothesis is physically realizable. In this case, the 
toroidal part would be equal to the said component of 
a“ and, as such, it would couple to the transverse elec¬ 
tromagnetic held. Nevertheless, in this hypothetical case 
the statement that the toroidal part couples to the held 
means nothing else than that a“ couples to the held, and 
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it would not amount to separate coupling of the electric 
and toroidal parts. 

The same conclusion applies to any value of j > 1. 
The root cause is the breaking up of the spherical Bessel 
functions in the split US- Their lowest order term, 
proportional to (fcr)- ,_1 , goes to the electric part and the 
rest of the terms go to the toroidal part. When unsplit, 
the spherical Bessel functions in exact integral expres¬ 
sions like Eq. © or Eq. (El completely reject the 
|p | ^ w/c components of J u (r). The Fourier transforms 
of spherical Bessel functions correspond to momentum 
space radial deltas <5(|p| — w/c) (see end of Sec. Ill in 
Ref. ED). After the split, the term (fcr)- 7 1 by itself does 
not provide such rejection. This allows |p| ^ ui/c compo¬ 
nents to “leak into” the electric part, and to be forcefully 
present in the toroidal part as well since the sum of the 
two parts does not contain any out of shell contribution. 

We have hence established that the electric and 
toroidal parts cannot be separately determined by mea¬ 
suring the electromagnetic fields produced by the source 
outside its region, or by measuring its coupling to ex¬ 
ternal electromagnetic waves. Our result implies that, 
as opposed to electric and magnetic multipoles whose 
distinct character is given by their parity, there is no 
electromagnetic field of pure toroidal character. It also 
implies that toroidal multipoles cannot be independently 
excited, i.e. there is no independent toroidal coupling. 
This precludes the existence of selection rules of toroidal 
character. The two different parities (electric and mag¬ 
netic), and their linear combinations are necessary and 
sufficient to characterize electromagnetic transitions (see 
e.g. Chap. 9.8 p. 436 in Ref. EE and Chap. XII 2. B 
in Ref. I26T ). The lack of longitudinal coupling together 
with the binary character of the parity operator leave no 
room for a third option. 

While the a“ TO have an independent physical meaning, 
the individual electric and toroidal parts do not have it. 
The toroidal multipoles cannot be considered an inde¬ 
pendent family. The same can be said about the isolated 
electric part, or about any split of the terms in the Tay¬ 
lor series of the spherical Bessel functions in Eqs. (E9) 
or Eq. ©■ Any such split is bound to introduce non- 
radiative out of shell components in each of its parts. 

Since the toroidal multipoles are higher order terms in 
the approximation of physically meaningful quantities, 
they can be used to improve analytical models that use 
only the lowest order terms. The inclusion of the next to 
leading order term improves the accuracy and prediction 
ability of such models, as can be seen in Refs. 0, he m 
[TIM In Ref. E| and in the context of non-radiating 
source configurations, the inclusion of the second order 
term explains a zero in the induced , which cannot be 
exp lained by the first order term alone (see Fig. 2 in Ref. 
[20j. The improvement of the models will be maximized 
by using expressions like Eq. ©, Eq. (20) in Ref. ED . 
and Eq. (IA2I) . These expressions inherently contain all 
correction orders because they are exact for any size of 
the source. 


C. Splits without out of shell components 

We have seen that the out of shell non-radiative com¬ 
ponents appear when splitting a“ m by separating terms 
of different orders in the size of the source. Nevertheless, 
we observe that the a“ m in Eqs. (IA2I) [a“ m in Eq. ©] 
have two formally distinct contributions, and that each 
of them contains a spherical Bessel function which guar¬ 
antees the restriction to |p| = uj/c, thereby avoiding the 
intrusion of the out of shell components. However, we 
will now show that, due to the absence of longitudinal 
radiation, the two on-shell contributions cannot be sep¬ 
arately determined. We consider another orthonormal 
basis for vector functions on the momentum shell: The 
vector spherical harmonics Yy; irn (p). They are related 
to the {Xj- m (p), Z jm (p), W im (p)} as 

Xy m (p) = Yjj im (p), 

Zjm(p) = — y/^fYjJ_l, TO (p) — 2j+l Yj,j'-|-l ) m(p), 

W im (p) = Y jJ+1 , m (p).(7) 

The first and third lines of Eq. m can be found in 
Eq. (16.91) and Exercise 16.4.4 of Ref. [El, respec¬ 
tively. The second line can be obtained from the equa¬ 
tion "Zijm^/jij + 1) =ipx ~LYj m . The result is reached 
by using the expression of the orbital angular momentum 
operator in momentum space L = «V p x p, the formula 

px (VpXp) = —(P'Vp)p, and the gradient formula 

Eq. (16.94) in Ref. 0 

Figure [2] illustrates Eq. ©. Tabic [H] contains informa¬ 
tion about the Yj ! ; jm (p). Notably, we deduce from Eq. 
© and see in Fig. [2] that, while Yjj im (p) is transverse, 
Yjj_i im (p) and Yjj+i, m (p) are of mixed character, i.e., 
neither transverse nor longitudinal. 


J 2 

Jz 

n 

L 2 

Polarization 

Yjj,m(p) j(j + 1) 

m 

(-i )* +1 

j{j + 1) 

X 

(p) j(j + 1) 

m 

(-i)' 

(j - Y)j 

mixed 

Yjj + l,m(p) j(j + 1) 

m 

(-ir 

{j + l)(j + 2) 

mixed 


TABLE II. Vector spherical harmonics: Polarization charac¬ 
ter and eigenvalues of J 2 , J z , n, and L 2 . L 2 is the orbital 
angular momentum squared. Both j > 0 and m = — j ... j 
are integers. For j = 0, Yo,i,o = — p/(4-7r), Yo,o,o = 0, and 
Yo,-i,o is not defined. 

We can expand J w (p) as: 

JUp) = (8) 

+ T j,j — l,rn^j,j-l,rn(p) + T j,j+l,rn^j,j+l,m{p)- 

j,m 
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It follows from Eq. 0 that: 


■ 

jm 


.7 + 1 

2j+l 


2j+l 

We invert Eq. 0 to obtain: 


2j + l 


1+1 -u 
2j+l 'j,j+l,m- 


(9) 


T u . =b“ 


j,j- 1,™ 




1+1 


2j+l “jm 


a + 

U/ nm. I 


a,™ — 


2j+l 




J + l ^ 

2 j+l L 'jm‘ 


( 10 ) 


The current dependent integrals in r space that define 
m contain only the spherical Bessel function of order 
l. This can be shown using Eq. m and the expressions 
for {a" m ,6“ m ,c“ m } in Eq. (1A2[) . For example, in the 
dipole case of Eq. 0, the integrand that contains jo{kr) 
is only due to r“ 0 m , and the integrand that contains 
j 2 (kr) is only due to Tf 2 ,m- Both terms have zero out of 
shell components. 

Let us assume that we want to determine r“j_ ± m and 
T j’,j+i,m by measuring the electromagnetic radiation from 
the source at points r such that |r| > R. As is clear 
from Eq. m, and can be visually appreciated in Fig. 
[2] the T^j ±l rn contain both transverse and longitudinal 
components. Their measurement involves the detection 
of both transverse and longitudinal fields generated by 
the current distribution. But, as previously discussed the 
longitudinal field due to the current distribution outside 
the source region is exactly canceled by that due to the 
charge distribution. Without access to the longitudinal 
degrees of freedom of J w (p), it is impossible to separate 
the two Tj’j .|_i m contributions to af m . 

We see that there are splits of af m that, while both 
parts remain free of out of shell components, they con¬ 
tain longitudinal components which render them non- 
separablc. This is the case for a definition of the toroidal 
multipoles that has been recently given in Box 2 of Ref. 
Il9l . In App. [0 we first show that this definition is dif¬ 
ferent from the original definition in Refs. @ and 0. In 
the split of the a“ TO proposed in Box 2 of Ref. m the in¬ 
tegrands contain entire spherical Bessel functions which 
prevent the appearance of out of shell components. We 
then show that each of the two parts contains longitu¬ 
dinal components that cannot be detected from outside 
the source region. 


IV. DISCUSSION AND CONCLUSION 

To finalize, we discuss some statements that are often 
found in the toroidal literature. 

It is often stated that the toroidal multipole family 
is a third family, independent of the electric and mag¬ 
netic ones i a El EH • Originally, the authors of Ref. 


@ suggested in Sec. 4 to replace the {af m , &“ m , c“ m } 
by {tj m , 6" m , dj m }, where are the coefficients of the 
toroidal part of af m . We have shown that the i“ TO cannot 
be separately determined, and can hence not constitute 
a third independent multipolar family. 

It is also often stated that three multipolar families, 
including the toroidal one, are needed in order to expand 
a general charge-current distribution, while the radiation 
that they emit can be described with only two multipolar 
families seiei. Starting from Eq. 0. our analysis 
makes clear that, indeed, three multipolar families are 
needed to describe the source, and only two are needed to 
describe the fields radiated by it. It also makes clear that 
the third family is constituted by the longitudinal multi¬ 
poles, and not the toroidal multipoles. The fact that the 
longitudinal field is zero outside the sources explains why 
the number of necessary multipolar families decreases by 
one when going from the sources to the fields. 

Finally, it is also stated in Refs. HI and El that the 
toroidal multipoles can be separated using spectroscopic 
techniques. The separation is to be enabled by exploiting 
the extra frequency dependent weighting of the toroidal 
multipoles [ k 2 = (w/c) 2 ] with respect to the electric mul- 
tipoles. It is straightforward to prove that, should this 
be actually possible, there would exist a doubly infinite 
number of independent multipolar families, and not just 
a third one. This unphysical outcome shows that the 
conjectured separability is not possible. In order to show 
this, let us assume that it is possible. The steps going 
from Eqs. 0 to 0 reveal that the extra k 2 scaling 
comes from different terms in the series expansion of the 
spherical Bessel functions. We stopped at the second low¬ 
est order, but each of the infinite number of terms in the 
complete expansion has an extra k 2 weight with respect 
to its predecessor. Additionally, one can make a simi¬ 
lar expansion of the magnetic multipoles (see Sec. 4 
in Ref. l2lh . If the separate measurement of the terms 
with additional k 2 factors where possible, as conjectured 
in Box 2 of Ref. EH there would be not just an ex¬ 
tra independent multipolar family, but a doubly infinite 
number of them. The physical reality is that each of this 
doubly infinite number of formally obtained terms is not 
an independent degree of freedom. The multipolar elec¬ 
tromagnetic transitions that determine the spectroscopic 
response of a system of charges and currents are com¬ 
pletely characterized by using the two possible parities, 
electric and magnetic. The spectroscopic measurements 
of {a" m ,6" TO } already contain all the available informa¬ 
tion. On the other hand, should the analysis of the spec¬ 
troscopic data make use of models limited to the lowest 
order, the inclusion of the second lowest order terms can 
only improve it, and the use of exact expressions like Eq. 
0, Eq. (1A2[) . and Eq. (20) in Ref. Hi] will optimize it. 

In conclusion: Our analysis proves that the dynamic 
toroidal multipoles do not have an independent physical 
meaning with respect to their interaction with electro¬ 
magnetic waves. Their formal meaning is clear, however: 
They are higher order terms of an expansion of the trans- 
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Appendix A: Contributions to 6“ m , o" m and c" m where ji(-) are the spherical Bessel functions of the first 

kind. 

Equation (14) in (2l} is an exact expression for the 
{a“ m , Vjmi c j m } coefficients in terms of integrals in both 
momentum and coordinate space. With Q jm (p) stand¬ 
ing for any of the {Xj m (P) ? Z jm (p),W jm (p)} and 
for any of the corresponding {a“ TO , b“ m , c" m }, Eq. (14) in 
HU reads 




As shown in [2l|, App. B], only terms with l = j con¬ 
tribute to the while the a^ m and c^ m get contribu- 

dp Q] m (p)n™(P) [ d 3 v J w (r)^(f)jKfcr), tions from both 1 = 3 - 1 and ~ l = 3 + L Explicitly: 

L 1 


(Al) 


— ™=j j- 

^ b jm = (~ i y E J d P X jm(P) Y Sm(p) J d 3 r J u (r)Y* m (r)jj(kr), 


m=—j 


-a“ = 

2 om 


m—j — 1 


i~i) j 1 X] [dp Zj jm (p)Y j _ lm (p) f d 3 r J ul (r)Y*_ lm (r)jj_ 1 (kr)+ 
m=j+l 

( -i) j+1 E J d P Z Um i+ ^(p) J d 3 r J u (r) 17 + 1 ?B (r)j i+ 1 (fcr), 


m=-0'+l) ' 


-c- = 
2 


m=j-l 


r r 

H ) j_1 E J dpW] rn {p)Y j _ 1Wi (p) J d 3 r 3 UJ (r)Y*_ lm (r)j j _ 1 (kr)+ 

1) 

™=J + 1 j. „ 

H ) i+1 E J dp W] m (p)Y j+lm (p) j d 3 v 3 u (r)Y; +lm {T)j j+1 (kr). 


m=-(j+ 1 ) ' 


(A2) 


Appendix B: Cancellation of longitudinal fields 
outside the source 

We show that the longitudinal field with |p| = u>/c is 
zero outside the source region. 

We consider spatially confined monochromatic electric 
charge and current density distributions p u { r) and J w (r) 
embedded in an isotropic and homogeneous medium with 
constant permittivity e and permeability p. We assume 
them to be confined in space. In the Lorenz gauge, the 
scalar and vector potentials meet the following inhomo¬ 
geneous wave equations: 

(V 2 + (ca/c) 2 ) r)= P ^ T \ (B1) 

(V 2 + (w/c) 2 ) A w (r) = -/iJ w (r), 

where c = 1/^/e/I. Outside the source region, it can be 
shown that the spatial Fourier transforms of (f> ul ( r) and 


A w (r) are non-zero only for |p| = w/c (see [U App. A]). 
With the help of [ 2 ^, Eq. (3.8)], they are readily seen to 
be proportional to the Fourier components of the sources 
in the same spherical shell domain: 

k P) = A(P) = (B2) 

In coordinate space, the electric held as a function of 
the potentials is 

E w (r) = iwA^(r) - V</> u (r), (B3) 

which, in momentum space (V —> ip) reads 

E u (p) = iwA(p) - ip</>(p). (B4) 

The longitudinal electric Held is hence 

P E w (p) = p (iwA(p) - ip&ipf) . (B5) 







Using Eq. (IB2I) and that p is restricted to p = we 
can write Eq. (IB5I) as 
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47 rk 


Esca(l') — _ ^ ' [Q jm(r) + Mj m ^jm(^) “I - 5 




u; ^ ^ ^, c<j ^ 

-p E^(p) = -p 

c c 


■ M*Mp) 

• pUP)\ 


c f 

^ 4t r 

iP 47TC ) 

(B6) 

sMTT)] 


lU1>L (^P • j«(p) - w/i;(p 


47T 


Mjm — 


Witi +!) 


d 3 r Pul (r)Y* m (r) — [rjj{kr )] , 

J d 3 r [V-(r x J w (r))]y/ m (f)jj(fcr), 


The term inside the brackets in Eq. (IB6I) is equal to 
zero because of the continuity equation in momentum 
space 


V • J w (r) = ico Pul ( r) ip • J w (p) = iuj Pul {p) (B7) 

particularized at |p| = w/c. 

i-P ■ -L(p) = *wp w (p). (B8) 

c 


We conclude that, outside the source region, the lon¬ 
gitudinal field with |p| = ui/c produced by the current 
density exactly cancels the one produced by the charge 
density. Note that the result is gauge independent since 
it is a statement about the electric field. 


Tjm = yfj(j + 1 ) / d3r ^ ‘ Y J*rn(*)jj( kr )> 

(Cl) 

where E sca (r) is the field produced by the sources, and 
(r) and $ jrn (r) are the multipolar fields of electric 
and magnetic parity, respectively. The Qj m are said to be 
charge excitations yielding electric multipoles, the Mj m 
transverse (w.r.t r) current excitations yielding magnetic 
multipoles, and the T lrn radial current excitations yield¬ 
ing toroidal multipoles. 

We first note that Q ]rn and are both multiplying 
the same multipolar field of electric parity VH (r) . This 
implies that, together, Qj m and Tj m must completely de¬ 
termine the multipolar coefficients of electric parity. This 
can be readily checked by setting the magnetic currents 
to zero in Jackson’s [27l . Eq. 9.167] expression for the 
exact multipoles of electric parity a^(j, m), namely: 


aE(j,m) 


k 2 

Witt + 1) 


Appendix C: An alternative definition of the 
toroidal multipoles 


In this appendix we first show that the definition of 
toroidal multipoles recently given in 19, Box. 2] is dif¬ 
ferent from the original definition in |2j, [3j. Both defini¬ 
tions involve the split of the multipoles of electric parity 
into two parts, but the splits are different in [l9|, Box. 
2] and M- We demonstrate this by showing that the 
well known expression of the toroidal dipole in the limit 
of small source 3|, Eq. 2.11] cannot be recovered from 
the definition of |19j, Box. 2]. Some terms are missing. 
We then also show that the missing terms are contained 
in a coefficient which is explicitly excluded from the def¬ 
inition of toroidal multipoles in [191 Box. 2]. The deriva¬ 
tions in this appendix recover one of the results from the 
main text: The toroidal dipole is just the next to lead¬ 
ing order term in the small source expansion of the exact 
electric dipole. The difference is that here the result is 
obtained directly in coordinate (r) space, while Eq. © 
from the main text was obtained in [ 211 ] by first going to 
momentum (p) space, and then going back to r space. 
Finally, we show that, in this alternative definition, both 
parts contain longitudinal terms, which render them non- 
separable. 

We start by examining the definitions in [l9], Box. 2] 


J d 3 v Y* m ( f) |cp w (r) ^ [ rjj{kr )] +ik[ r • J w (r)] ^-(fcr) j . 

(C2) 

It is clear from Eq. (IC1I) that Qj m corresponds to the 
first term of the sum in Eq. (IG2|) . and iTj m to the sec¬ 
ond term. The sum Qj m + iTj m determines a,E(j, to). 
The definition of Tj m in Eq. (IC1I) involves a split of the 
a e (j, tti) into two parts. We now show that it is a differ¬ 
ent split from the one in the original definition of toroidal 
multipoles 111 - 
Let us use iTi m 
nh r 

iTlm = 71 J d3r ^ ‘ Y * m ^i(kr), (C3) 

to attempt to recover the toroidal dipole in the small 
source approximation a Eq. 2.11], 

*1 =/ d 3 r ^{[r-J w (r)]r-2r 2 J aj (r)}, (C4) 

and which follows from the original definition (II Eq. 
39], 0, Eq. 2.11]). 

We start by arranging the three components corre¬ 
sponding tom = {l,0,—1} into a vector, 


• T s p h 


T 11 

T 10 

Ti_i 


ik 

71 


d 3 r [r • J w (r)] 


2 i: 

Y, 


10 


I i_i 


3i(kr), 

(C5) 
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and consider the correspondence between Y* m (r) and r 
in the spherical basis 


r 

h 

„ /tt 


„ [tt 

v* 

Ml 

i r i 

1 

_ 1 

= 2 J- 
V 3 

_1 

= 2,/- 
V 3 

1 

A 2 *1 


Equation (1C5I) is a vector in the spherical vector ba¬ 
sis except that, when both r and J u (r) are expressed in 
the spherical vector basis, the term [r • J u (r)] should be 
written [/J w (r)], where r' denotes the hermitian con¬ 
jugate of the position vector r. This is due to the fact 
that, in the spherical vector basis, r has complex com¬ 
ponents. In the Cartesian basis used in [l^, Box 2] and 
Eq. (EH), r is real valued and the dot product [ivJ w (r)] 
can be written as [r • 3 U (r)] . In this appendix we will use 
the Cartesian basis from now on. We can change Tj ph 
to the Cartesian basis by multiplying T® ph itself, and the 
other vectors involved in the expression with the change 
of basis matrix 


Ojx 


r^i 0 -4=1 

V2 U V2 


a 1 

dy 

= 

0 =±- 

y/2 y/2 


a 0 

a z_ 


O 

T - 1 

O 


a~ 1 


(C7) 


After using Eq. (IC6I) and changing the basis, the 
Cartesian expression reads: 

*Ti = ^ J d 3 r [ r ' -Mr)] r^^j^kr). (C8) 

We now use ji(kr) rs kr/ 3 in the limit of small fcr, to 
approximate Eq. (IC8I) by 


*Tr 



d 3 r [r • 


J 4 r )] 


„ fz 1 kr 

r V 7T 2 ”3" 


ik 2 /If 


3 [r-J„(r)]r 

6 


(C9) 


Equation (IC9I) cannot reproduce Eq. (ICdl) because the 
r 2 3 u (r) terms present in Eq. (IC4I) are missing in Eq. 
(EH). The definition of Tj m in [lf| Box. 2] is hence 
differen from the original definitions of the toroidal mul- 
tipoles in 0S- 

We now show that the missing terms are contained in 
Q lm : They are the terms of second lowest order in the 
small source approximation of Qi m . This reproduces our 
results from the main text. 

We start with 


Q 1 m — 





(CIO) 


take the same steps as before regarding vectors and basis, 


Qi — 


a/2 . 


J d 3 r p u ( r )r^~-!j- [rj^kr)], (Cll) 


and use the continuity equation p u (r) = V • (r) / (ike ), 

to obtain: 


(C12) 

We now consider the term in the shaded box, which, 
using the derivative of spherical Bessel functional can be 
written as 

[rj 1 (kr)\ =j 1 (kr) + jy \jo(kr) - 2j 2 (fcr)]| • (C13) 

We now take terms up to order (kr) 3 in the small argu¬ 
ment approximation of the spherical Bessel functions of 
Eq. (ICTall : 


jo (kr) R 

s 1 v 

ji(fcr) r 

kr ( 

w( 

32 (kr) r 

(kr) 2 

15 


(kr) 2 


1 - 


(kr) 2 

10 


(C14) 


to obtain the approximate expression 
d 


d [rji(kr)\ « 
kr / (kr) 2 

y l 1 ltT 


kr 

y 


\ _ (kr f _ (kr) 2 

6 15 


(C15) 


We now plug Eq. (IC15I) into the right hand side of Eq. 

EH 


Qi 


-i 3 1 


kV2\ 7T 2 

- i [31 f 


d 3 r [V-J„(r)]r^ + 


kV2\ Ti- 2 


d 3 r [V • J u (r)] 


■ (^’ r ) 


1 1 2 
10~6"15 
(C16) 


and reduce it to 

— l 


Qi 


V&7V 


= J d 3 r [V • 3 U (r)] r+ 




[V ■ J/r)]r 2 r 
15 


(C17) 


The first line of Eq. (IC17I) can be readily brought^ to 
the familiar form of the approximate electric dipole for 
small sources: 


1 ^M x ) = 27TT dh- lW - (l + k)ji+i(x)) 

2 See the steps in l27t Eq. 9.14] 







































12 


i 


\f&TT 


/ 


d 3 r J w (r). 


A = J 0J (r)r 2 r k , and choose V as a spherical volume en¬ 
closing the sources, so that 3uj(r)r 2 r k is zero on its sur¬ 
face boundary. Then, the right hand side of Eq. (IC21I) 
(C18) vanishes. 

Therefore, using Eq. (IC19D for k = {1, 2, 3} we obtain 


This lowest order term is precisely the one separated 
by Dubovik 0,0 in his split of the exact multipoles of 
electric parity between what is referred to as the “elec¬ 
tric” part, which is this term in [23, Eq. 9.170], and the 
toroidal part, which are the higher order terms. 

Incidentally, reversing the use of the continuity equa¬ 
tion in the first line of Eq. (IC17D recovers the first line 
in [13, Eq. 9.170]. 

We are now interested in the next to leading order 
terms of Qi, which are shaded in the second line of Eq. 
(1C17I) . We now show that, when summed to Eq. (IC9I) . 
the toroidal dipole is recovered. The key step is to use 
the divergence theorem in the integration by parts of Eq. 
(IC17I) . To that end we first consider the following steps, 
where n and k run over {1, 2, 3} and S n k is the Kronecker 
delta: 


V • (Jr 2 r fe ) = ^2 [(d n J n ) r 2 r k + J n dn(r 2 r k )\ = 

n 

y. [(d n Jn) r 2 r k + Jn(2r n r k + r 2 S nk )\ = ( C19 ) 

n 

(V • J) r 2 r k + 2 (r • J) r k + J k r 2 . 


Since J u (r) is, by assumption, bounded in space, the 
divergence theorem can be used to show that 


J d 3 r V • [J w (r)r 2 rfc] = 0 for all k. 


(C20) 


/ d 3 r [V • J u (r)] r 2 r = 

J , , (C22) 

~ 2 Jd 3 r [r-J w (r)]r J d 3 rr 2 J u (r), 

whose left hand side appears in Eq. (IC17D . After substi¬ 
tuting Eq. (IC22I) into Eq. (IC17D we get: 

*71 \/f / d3r Y^{“ 2 [ r ' J ^( r )] r — ^ 2j ^(r)} ■ (C23) 
Let us recall that Eq. (IC'231) are the next to leading 
order terms in the small source expansion of Qi m - After 
summing them to Eq. (ICiil) 


^ 'll I dh ' ("I + 0 |r ■ J " (r)l r ■ ir 2j " (r) = 

J d 3 r { [r • J w (r)] r — 2r 2 (r)} , 

(C24) 

we recover the exact form of the toroidal dipole in Eq. 

(EH). 

This derivation also recovers a result from the main 
text: The toroidal multipoles are just higher order terms 
in the small source approximation of the exact multipoles 
of electric parity. 

To conclude, we note that, according to Eq. (1C12I) . 
which contains V- J w (r), Qi depends on the longitudinal 
degrees of freedom of J w (r). This follows from the r 
domain to p domain correspondence 

V • J u (r)-Mp ■ J w (p), (C25) 


The divergence theorem [13, p. 36] 


j d 3 rV-A = | A-rdS 


(C21) 


relates the integral of the divergence of any well-behaved 
vector field A over a volume V to the flux through 
the surface boundary of V. For our purposes we set 


which makes clear that the coefficients multiplying 
the longitudinal multipoles contribute to Qi [see Eqs. © 
and ©]. 

Since the sum Q^ m + iT^ m is equivalent to the fully 
transverse a^ m , it then follows that the longitudinal de¬ 
pendence introduced by Qi must be canceled by iTi. 
The split introduces longitudinal components in both 
parts, which renders them non-separable according to the 
discussion in the main text (see Sec. UVD . 
























